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THREE-DIMENSIONAL VISCOUS DESIGN METHODOLOGY 
FOR ADVANCED TECHNOLOGY AIRCRAFT SUPERSONIC INLET SYSTEMS 


Bernhard H. Anderson 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


SUMMARY 

The present report discusses a broad program to develop advanced, re- 
liable, and user oriented three-dimensional viscous design techniques for 
supersonic Inlet systems, and encourage their transfer Into the general user 
community. This program Is fourfold In nature, namely (1) develop effective 
methods of computing three-dimensional flows within a zonal modeling metho- 
dology, (2) ensure reasonable agreement between said analysis and selective 
sets of benchmark validation data, (3) develop "user" orientation Into said 
analysis, and (4) explore and develop advanced numerical methodology. 


INTRODUCTION 

Calculation of the extremely complex multidimensional viscous flow fields 
found In supersonic Inlet systems proposed for advanced technology aircraft 
presents difficult but necessary challenges. Although one method of determin- 
ing a suitable design system would focus upon an extensive testing program, 
the costs of hardware fabrication and Individual tests are such that an alter- 
nate, more efficient design procedure Is required. One method which Is cur- 
rently being pursued at NASA Lewis Research Center focuses upon the development 
of a series of computer codes to predict the aerodynamics of supersonic Inlet 
systems In general geometries operating at realistic flow conditions. When 
developed further, these analyses could be used to (1) examine basic flow 
mechanisms governing the flow field under study, (2) give guidance to new and 
more efficient component designs, (3) predict the system operation at various 
flight conditions, and (4) reduce the experimental test matrix. All these 
Items would be of significant value to the design engineer. 

The broad objective of the program at NASA Lewis Research Center Is to 
develop advanced, reliable, and user oriented computer design techniques for 
advanced supersonic Inlet systems and encourage their transfer Into the general 
user community. 

(1) Develop effective methods of computing the three-dimensional viscous 
flow fields within supersonic Inlet systems using Spatial Marching and Navler 
Stokes techniques within an overall zonal calculatlonal methodology. 

(2) Generate a significant benchmark experimental data base against which 
said analyses may be verified; determine the range of applicability and accur- 
acy of said analyses, and provide guidance for Improving these analyses by 
comparison with benchmark experimental data. 




(3) Develop “user orientation" Into these analyses which Includes graph- 
ical output representation, case "running" protocol for flows of Interest, and 
a formalized structure to use these three-dimensional codes. 

(4) Explore the development of Improved three-dimensional grid genera- 
tion, numerical methodology, and computer technology utilization relative to 
supersonic inlet analyses. 


ZONAL ANALYSIS METHODOLOGY 

Generation of analyses for predicting the viscous flow In advanced super- 
sonlc Inlet systems represents a difficult but Important problem. The flow 
field In physically realistic situations Is three-dimensional, viscous (and 
usually turbulent), and contains strong shock waves. One method of developing 
an analysis for these flow fields would focus upon the solution of the three- 
dimensional Navler Stokes equations. Although these equations could be used 
In principle for the entire flow field of Interest, the number of grid points 
required to obtain adequate flow resolution Is prohibitive and makes this 
approach Impractical even for the next generation of computers. 

A suitable alternative which has been pursued by the Aerodynamic Analysis 
Section at the NASA Lewis Research Center Is the zonal embedding methodology. 
Under this philosophy the overall flow field Is divided Into distinct parts 
and each part analyzed via the appropriate set of equations. Interaction 
between sections may be considered as required. The use of zonal embedding 
methods Is advantageous In several ways, namely (1) the region over which the 
complex equations are used Is minimized, (2) the simpler sets of equations 
have faster algorithms which result In Improved computational efficiency and 
(3) the computer storage requirements for the overall procedure are not as 
great. All of these advantages lead to more cost effective computer analyses. 

An example of the zonal embedding methodology is Illustrated In figure 1 
which Is developed for three-dimensional supersonic Inlet systems. The flow 
In the supersonic portion of the Inlet will be calculated using the three- 
dimensional PEPSIS general flow field solver. Downstream of the Inlet throat 
In the subsonic portion of the Inlet diffuser, the flow will be computed using 
the PEPSIG general flow solver. Interfacing these two regions Is the Inlet 
throat. Within this region, there Is supersonic flow upstream of the normal 
shock wave and subsonic flow on the downstream side. In this region of the 
Inlet, the flow field will be analyzed using the three-dimensional generalized 
MINT computer code. The computations In this region would be numerically 
Interfaced with the upstream PEPSIS calculation and the PEPSIG downstream so- 
lution. A similar zonal methodology was developed for the analysis of forced 
mixer nozzles and nonaxlsymmetrlc nozzle systems. 

In 96^ era l* the flow fields of Interest can be analyzed via one or more 
of the following approaches: (1) supersonic three-dimensional viscous spatial 
marching solvers, (2) subsonic three-dimensional viscous spatial marching 
solvers, and (3) three-dimensional Invlscld/vlscous elliptic solvers. General- 
ized analyses for each of these approaches have been developed and validated 
against both laminar and turbulent benchmark experimental data. These General- 
ized Three-Dimensional Flow Solvers are called the PEPSI/MINT series of com- 
puter codes and are described In the appendix. The PEPSIS computer code 
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(refs. 1 to 5) was developed as a generalized three-dimensional supersonic 
viscous marching solver; the MINT code (refs. 6 to 16) as a general Invlscld/ 
viscous elliptic flow solver; and the PEPSIG code (refs. 17 to 21) as a gener- 
alized subsonic three-dimensional viscous marching flow solver for diffuser 
ducts. The PEPSI/MINT series of viscous flow solvers are "general" because 
the geometry, computational mesh along with all the necessary metric Informa- 
tion, and the boundary conditions are external to the solution generator, l.e., 
this Information must be supplied as Input data. 

Each of these three-dimensional flow field solvers has proven to be a 
successful technique for solving the set of equations at relevant flow condi- 
tions In relevant geometric configurations. They represent efficient and 
accurate methods for analyzing advanced supersonic Inlet systems. In particu- 
lar, the PEPSI/MINT series of general three-dimensional flow solvers was writ- 
ten In modular form so that changes In Items such as boundary conditions, 
turbulence models, and geometry as well as additional terms representing new 
physical phenomena can be made In a straight forward manner. Changing one com- 
ponent, such as geometry or form of the governing equations, would not require 
changes In other parts of the code, such as the subroutines which solve the 
sets of equations. The PEPSI/MINT series of generalized flow solvers has been 
designed to provide the maximum flexibility to anticipate future needs and 
respond to those needs. Analysis of new configurations will thus focus upon 
the development of the geometry and mesh and In developing smoothness criteria 
for the particular flow algorithm and fluid dynamic characterlsts of the flow 
Itself. Since these three-dimensional flow field codes have been extensively 
verified for a selected set of relevant benchmark experimental data (refs. 22 
to 53) they represent the general set of spatial marching and elliptic solvers 
which forms the basis for ongoing and future computational efforts. 


BENCHMARK VALIDATION AND IMPROVEMENT 

The goal of providing design analyses tools for supersonic Inlet system 
technology advancement necessitates a two tier experimental validation process 
to properly evaluate the ability of these three-dimensional solvers to predict 
the relevant flow physics. First, It Is necessary to show good predictions of 
flow parameters of Interest to the user community on "generic" configurations 
which would be typical of real world designs. Such testing should be as 
realistic as possible. Confidence In the advanced computational capability 
would Increase by leaps and bounds In the user community by this type of 
demonstration . 

A second tier program appears to be called for and here essential ele- 
ments of these advanced computer design tools would be validated in very simple 
environments that highlight one or more basic flow mechanisms. The require- 
ments are for overall simplicity and easy access to make the necessary very 
detailed measurements. Detil validataion of basic flow mechanisms is required 
because the generic component tests may or may not emphasize any particular 
basic flow mechanism and overall good agreement may be just that, an "overall" 
agreement. There may be an unfortunate cancelling effect In the basic flow 
mechanisms of generic experiments. For Instance, one could not analytically 
Investigate the factors affecting Inlet performance with any degree of confi- 
dence If the Individual shock wave turbulent boundary layer Interactions en- 
countered within the Inlet were not well predicted (fig. 2). Confidence In 
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making extrapolations using computational tools will Increase If the basic 
mechanisms, when Isolated, are predicted both qualitatively and quantitatively. 

The PEPSI/MINT series of general flow field solvers are large, complex, 
computer analyses designed to solve the three-dimensional flow field charac- 
teristic of advanced aircraft Inlet concepts. As a result, benchmark valida- 
tion and generic verification of these codes Is a very Important part of the 
program at NASA Lewis Research Center. A series of experimental data sets 
have thus been Identified from the literature or obtained through university 
contracts or grants as well as Inhouse experiments to be compared with the 
PEPSI/MINT codes. 


Three-Dimensional Supersonic Viscous Marching Flow Solver (PEPSIS) 

Following this philosophy on computer code validation, a series of generic 
Inlet configurations have been Identified from reports published and have been 
analyzed using the PEPSIS three-dimensional viscous flow solver. These super- 
sonic Inlets Include the Mach 2.5 mixed compression Inlet of Fukuda, Hlngst 
and Reshotko (ref. 22) the Mach 3.0 mixed compression two-dimensional Inlet of 
Anderson and Wong (ref. 23) the Mach 3.5 mixed axlsymmetrlc Inlet of Syberg 
and Hlckcox (ref. 24) and the Mach 7.5 hypersonic Inlet of Gnos and Watson 
(ref. 25). It was concluded In this series of generic inlet verification cal- 
culations that "overall" good agreement between experiment and analysis was 
achieved although a great deal of uncertainty exslsted In understanding the 
details of the bleed Interaction. This uncertainty centered around five main 
points; namely (1) the bleed distribution through the bleed region was unknown, 
(2) the size of the bleed holes was often large In comparison to the boundary 
layer thickness, (3) the bleed mass flow was often large In comparison to the 
boundary layer mass flow approaching the bleed region, (4) the possibility 
that recirculation exslsted In the bleed region, and (5) transpiration took 
place using discrete holes rather than continuous bleed. As a result, a major 
benchmark bleed program was Initiated at Lewis Research Center to obtain the 
necessary experimental data and physical understanding of the boundary layer 
bleed process. 

In addition to this set of calculations on generic supersonic Inlet con- 
figurations, a series of benchmark supersonic flow phenomenon has been Identi- 
fied and are presented In figure 3. These Interactions Include the oblique 
shock wave turbulent boundary layer Interaction of Rose (ref. 26) the three- 
dimensional glancing sidewall shock wave turbulent boundary layer Interaction 
of Oskam, V as, and Bogdonoff (refs. 27 to 29) and the cone-at-alpha interaction 
of Ralnblrd (ref. 30) . 

Two-dimensional oblique shock wave Interaction . - The first benchmark 
test case Is that of the oblique shock wave turbulent boundary layer Inter- 
action studied by Rose (ref. 26). The oblique shock was generated by a 10° 
cone placed In a blow down circular tunnel operating at Mach 3.88 and Reynolds 
number of 5.0E6 per ft. In this calculation, the entire tunnel was modeled 
using a computation mesh of 89X450 (40 500 nodal points) with approximately 30 
points In the boundary layer. On the Lewis CRAY I computer, this benchmark 
test case took 0.9 min of CPU time. A comparison of the calculated and exper- 
imental wall static distribution through the Interaction zone Is presented In 
figure 4 and demonstrates that the PEPSIS solver correctly predicts the two 
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stage measured static pressure generated by the oblique shock wave Interaction 
with the boundary layer. The abrupt Initial static pressure was generated by 
the Impingement of the Incident shock wave, and the more gradual downstream 
static pressure rise was caused by the formation of the reflected shock system. 
Figure 5 presents the comparison between the calculated and measured Mach num- 
ber profiles through the Interaction zone. The Incident shock appears as a 
discontinuity In the experimental data, while shock smearing occurs over sev- 
eral grid points. The PEPSIS analysis predicts well the development of the 
Mach number profiles Including the thinning of the boundary layer. A detailed 
discussion concerning the manner In which this calculation was performed 
appears In a paper by Benson and Anderson (ref. 37). 

Three-dimensional glancing sidewall Interaction . - A schematic diagram 
depicting the analytical and experimental test configuration used to study the 
glancing sidewall turbulent boundary layer Interaction Is shown In figure 6. 
This Interaction arises when the oblique shock wave formed by the wedge Inter- 
acts with the wall boundary layer. Because of the skewing of the Incident 
shock wave across the boundary layer on the tunnel wall, a strong transverse 
static pressure gradient Is established which generates strong cross flows. 

The nominal Mach number upsream of the wedge was 2.84, the tunnel total pres- 
sure was set at 88.9 psla and the tunnel free stream total temperature was 
445.0° Ranklne. Two computational mesh systems were used to study the glanc- 
ing sidewall boundary layer Interaction, a coarse grid composed of 40X40X100 
(160 000) nodal points and a medium grid with 40X60X90 (216 000) points. On 
the Lewis CRAY I high speed computer, these calculations were accomplished 
using 6.0 and 9.6 min of CPU time, respectively. Shown In figures 6(a) to (c) 
Is a comparison of the computed and measured yaw angle distribution through 
the wall boundary layer. The yaw angle Is defined as the ratio of the velocity 
In the YG-dlrectlon which Is parallel to the tunnel sidewall, divided by the 
velocity In the X-dlrectlon. This Is a particularly difficult parameter to 
calculate since It represents the ratio of two velocities which approach zero 
as the wall Is approached. It Is apparent that very good agreement was ob- 
tained using the PEPSIS analysis to model this Important Interaction. A detail 
study Into the effect of wall functions and mesh resolution was performed by 
Anderson and Benson (ref. 36) with the conclusion that Insufficient mesh reso- 
lution of the near wall region caused discrepancies to appear In the calcula- 
tion of the velocity field. 

A detail picture of the flow field In the tunnel sidewall region Is shown 
In figure 7. In the vicinity of the tunnel sidewall, a very strong vortex Is 
established which elongates and Increases In strength In the downstream direc- 
tion. This causes very low energy flow to accumulate In the corner region 
with resulting low wall shear stress In this region. This probably accounts 
for the very high heating rates measured by Oskam, Vas and Bogdonoff (refs. 27 
to 29). The Increased strength of the sidewall vortex Is also suggested by 
the experimental data results since the maximum yaw angle Increases In the 
downstream direction. The oveturnlng that occurs In the sidewall boundary 
layer results from an Imposition of the main stream static pressure gradient 
upon the low momentum near wall viscous flow. This overturning results In low 
energy fluid being drawn In towards the tunnel sidewall In the lower ramp 
region. In the outer ramp region, the flow rolls over to form the sidewall 
vortex. As the flow develops In the downstream direction, the sidewall vortex 
losses Its Identity and the surface velocity vectors tend to align themselves 
with the shock angle. 
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Figure 8 presents the static pressure signature calculated on the ramp 
and tunnel sidewall surfaces. The numerical formation of the shock wave In 
the region of wedge tip can clearly be observed as well as the shock wave for- 
mation on the tunnel sidewall surface. Although not explicitly shown In 
figure 8, an adverse static pressure gradient exslst In the near wall region 
of the wedge and tunnel sidewall surfaces. Of special Interest Is the fact 
that static pressure gradients are established on the ramp surface as a result 
of the corner flow boundary layer development. The static pressure gradients 
that are created on the ramp and tunnel sidewall surfaces fall Inside a coni- 
cal region of Influence. 

The global features of the glancing side Interaction that emerge from 
this study can be viewed as two separate flow regions, namely a large outer 
region of Invlscld high energy flow which follows the Invlscld stream lines, 
and a smaller Inner region near the wall composed of low energy fluid which 
migrates along the shock wave. The low energy flow will eventually accumulate 
In the corner region of the tunnel celling causing severe problems. This flow 
Is called the glancing sidewall/corner shock wave turbulent boundary layer 
Interaction and will be discussed later In this paper. 

Cone-at-alpha Interaction . - To understand the Importance of adequate 
mesh resolution In critical computational regions, a series of calculations 
were performed using the experimental configuration of Ralnblrd (ref. 30). 

This experiment consisted of a 12.5° half angle cone mounted In the supersonic 
wind tunnel at 15.6° angle of attack. The Mach number for this series of cal- 
culations was set at mach 4.25. The extreme angle of attack conditions were 
chosen because a recirculation vortex developes as a result of the Interaction 
of the windward and leaward boundary layers. 

Figure 9 demonstrates the Importance of radial and circumferential mesh 
resolution for this problem. The figure shows Mach number contours for the 
aft quadrant of the flow field for different meshes at a plane near the meas- 
uring station. All of the calculations required 220 streamwlse stations, with 
nearly 100 of them located near the cone tip to resolve the formation of the 
shock. Figure 9(a) shows the results of the first calculation with 49 radial 
points packed near the cone surface, and 19 circumferential points evenly 
spaced at every 10°. The results Indicate a thickening of the boundary layer 
on the leeward side, and give no Indication of a recirculation. Increasing the 
number of circumferential points to 37 evenly spaced at 5° while maintaining 
40 radial points gives the results of figure 9(b). This calculation gives an 
Indication of recirculation but with little detail. Because the recirculation 
Is caused by Interactions In the boundary layer, more resolution was required 
In the radial direction near the surface. Packing the 49 radial points closer 
to the surface gave the results shown In figure 9(c). While resolving details 
of the surface, resolution of the shock wave was sacrificed which resulted In 
violent post shock pressure oscillations. To resolve both the shock and the 
boundary layer, the radial mesh was Increased to 80 points. The results shown 
In figure 9(d) Indicate that the calculation was then able to resolve boundary 
layer details without oscillations near the shock wave. Figure 9(e) shows one 
further refinement In which the circumferential mesh was Increased to 50 points 
which were packed on both the windward and the leeward rays. The marked dif- 
ference In the results of figure 9(a) and figure 9(e) clearly demonstrates the 
need to resolve all of the small and large scale phenomena present In this 
problem. The 40X19X220 (167 000) mesh point calculation used 6.0 min of CPU 
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time on the Lewis CRAY I high speed computer, while the 80X50X220 (880 000) 
nodal point calculation used 31.7 min of CPU time. Figure 10 shows the com- 
puted recirculation vortex which formed on the leeward side of the cone flow 
field In terms of the secondary velocity vectors for the 80X50 transverse mesh 
system. Note the overshoot In the cl rcumf rentlal velocity In the region to 
the waterline side of the the recirculation vortex. 

In addition to these benchmark Interactions, a major experimental program 
Is underway at the NASA Lewis Research Center to obtain detail data on selected 
two and three-dimensional shock wave boundary layer interactions which focuses 
on the details of boundary layer bleed (fig. 11). Recent efforts to benchmark 
validate the PEPSIS flow solver have been directed towards understanding the 
details of these bleed Interactions and Improving the bleed model to accurately 
describe the flow physics of this Interaction. Data obtained by Hlngst and 
Tanjl (ref. 31) In the Lewis lxl Supersonic Wind Tunnel Is being used as the 
benchmark standard since detail boundary layer profile through the Interaction 
and bleed flow distribution through the holes were measured. Along with the 
shock wave boundary layer bleed Interaction study, the glancing sidewall/ 
corner shock wave Interaction, previously mentioned, Is also being examined 
experimentally In the 1X1 ft Supersonic Wind Tunnel. 

Glancing sidewall/corner Interaction . - Figures 12 and 13 present the 
analytical and experimental oil flow patterns resulting from the the Impinge- 
ment of the glancing sidewall flow field with the tunnel floor of the 1X1 ft 
Lewis Supersonic Wind Tunnel. The mesh system used to study the glancing 
sidewall/corner Interaction was composed of 80X40X120 (384 000) grid points. 
Calculations were performed on the Lewis CRAY I computer using 14.4 min of CPU 
time. Although the shock static pressure rise was not sufficient to separate 
the flow In the analysis on the center portion of the tunnel floor, there Is a 
corner region separation that was revealed In the calculations and can be seen 
In the experimental oil flow patterns. Because this Interaction Is typical of 
a class of problem corner Interaction that can be found In the advanced super- 
sonic Inlets concepts, It was deemed Important to obtain detailed experimental 
benchmark data In the Lewis 1X1 ft Supersonic Wind Tunnel. 

The PEPSIS three-dimensional viscous analysis is also being used to Iden- 
tify selected benchmark phenomenon that will become Important In the design of 
advanced aircraft Inlet concepts (fig. 14). These highly three-dimensional 
shock wave turbulent boundary layer Interactions center around the hypersonic 
corners Interaction of Crescl (ref. 32) the Intersecting wedge Interaction of 
West and Korkegl (ref. 33) and the skewed wedge Interaction of Settles, 

Perkins, and Bogdonoff (ref. 34) and the crossed side shock wave turbulent 
boundary layer Interaction for which no benchmark data is available. 

Hypersonic corner Interaction . - The primary interest in the experiment 
by Crescl (ref. 32) was directed towards the hypersonic, low density flow re- 
gime wherein the Invlscld shock layer and the viscous boundary layer produced 
under these conditions are of the same order of magnitude. In this case, the 
viscous and Invlscld effects are Interrelated and cannot be treated indepen- 
dently. This phenomena Is Important In that It represents a complex corner 
Interaction that Is typical of a large class of problems associated with ad- 
vanced supersonic Inlet concepts. The structure of the corner flow Interac- 
tion Is shown In figure 15. Qualitatively, the structure Is similar to that 
found by West and Korkegl (ref. 33) and consists of an oblique corner shock 
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resulting from the Intersection of the shock waves generated by the rapid 
boundary layer growth on the wall surfaces, and a set of two trlpple points 
formed by Intersection of the embedded shock wave system. This Interaction was 
studied by Buggeln, McDonald and Kin (ref. 5) using the PEPSIS forward march- 
ing three-dimensional flow solver, as part of the benchmark verification proc- 
ess. The model was mounted In a Mach 11.8 blowdown wind tunnel, and tests 
were conducted at free stream Reynolds varying between 0.15E6 and 0.5E6. The 
tunnel stagnation temperature varied between 1700° Ranklne and 1900° Ranklne 
which produced wall temperature ratios of 0.29 and 0.32. Calculations were 
performed on the hypersonic corner Interaction using a grid system composed of 
50X50X120 (300 000) nodal points. Although the calculations were performed on 
the Lewis IBM 370/3033, the equivalent computing time on the Lewis CRAY I 
computer would be 10.8 CPU min. Figures 16 and 17 present the calculated Mach 
number and total pressure field of the hypersonic corner Interaction. The 
primary features of an oblique corner shock wave generated by the Intersection 
of the the boundary layer shock wave from the wall surfaces and the set of 
trlpple points formed by the embedded shock system are easily seen In these 
figures. In addition, the secondary feature of a triangular region bounded by 
two slip surfaces and the corner shock wave was also revealed In the 
calculations . 

Lewis 40/60 mixed compression Inlet . - Figures 18 through 20 present a 
series of calculations on a Mach 2.5 mixed compression supersonic Inlet at 0° 
and 2° angle of attack. The purpose of the study was threefold, namely (1) to 
examine the problems associated with computing boundary layer bleed within an 
Inlet envlorment, (2) to study the effects of boundary layer bleed on Inlet 
characteristics, and (3) to anticipate flow problems associated with the three- 
dimensional shock boundary layer Interactions of this Inlet at angle of attack 
and recommend solutions. A detailed experimental bleed study was performed on 
this Inlet by Fukuda, Hlngst and Reshotko (ref. 22) In the 10X10 Supersonic 
Wind Tunnel at Lewis Research Center. This axlsymmetrlc mixed compression 
Inlet was designed for a Mach number of 2.5 with 40 percent external area con- 
traction and 60 percent Internal contraction. The external compression was 
accomplished with a 12.5° half angle cone and the Internal cowl lip angle was 
0.0°. Figure 18 presents the results of a two-dimensional calculation using 
wall functions to establish a baseline or reference case. The two-dimensional 
computational mesh was composed of 89 radial points and 360 forward marching 
steps for a total of 32 040 nodal points. On the Lewis CRAY I high speed com- 
puter, this calculation used 0.7 CPU min. Shown In figure 18 are the local 
Mach number profiles computed within the Inlet geometry. The cone shock as 
well as the reflected shock system within the Inlet are clearly vlsable. The 
cowling produced an Internal shock wave system which reflected first off the 
centerbody and then off cowling reaching the throat region for a second center- 
body Interaction. The second centerbody reflected shock was the terminal shock 
when the model was operated In the critical mode. 

As an Initial start to examine the three-dimensional shock wave turbu- 
lent boundary layer Interactions that can arise within axlsymmetrlc Inlets at 
angle of attack, two separate cases were considered using the PEPSIS three- 
dimensional flow solver. In the first case, the 40-60 Inlet was operating at 
the design Mach 2.5 condition and 2.0° angle of attack, and In the second case, 
the 40-60 Inlet was operating at Mach 3.0 condition, also at 2.0° angle of 
attack. In each case, the centerbody was In the design position for zero angle 
of attack operation. The design Mach number calculation revealed a large 
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region of subsonic flow on the leaward side of the first cowl reflection. 

This Is In subtantlal agreement with limited experimental data at this condi- 
tion. Since little Information could be obtained on the three-dimensional 
character of the shock wave turbulent boundary layer Interactions within this 
Inlet for the first case, the overspeed case was run. 

The Mach 3.0 angle of attack case used a computational mesh composed of 
49X19 nodal points In the transverse plane and 480 forward marching step for a 
total of 446 880 grid points. The computing time for this case was 16.9 CPU 
min on the Lewis CRAY I high speed computer. Figures 19 and 20 present the 
results of this calculation. For the Inlet operating at the design Mach number 
and at 0° angle of attack, the cone shock Is located just outside the cowl 
lip, thus spilling about 0.5 percent over the cowl. In the Mach 3.0 overspeed 
case, the cone shock Is forced Inside the cowl lip. At angle of attack, this 
shock system translates forward on the leeward side and rearward on the wind- 
ward side, thus creating a small expansion region on the cowl lip. Figure 19 
shows surface static pressure distribution on the cowl and cone. The shock 
waves are clearly shown as light bands which are Inclined relative to the Inlet 
axis. On the cone, the shock system moved forward on the leeward side and aft 
on the windward side, while on the cowl, the first rather weak shock from the 
cone tip and the stronger shock from the cone reflection are clearly evident. 

Also evident on the Inlet surfaces are pressure contours which are normal to 

the Inlet axis, particularly near the throat section on the cone and along the 
cowl between the shocks. These features are the results of compressions and 

expansions present In the Inlet due to surface contouring. For an other per- 

spective of this flow field, static pressure and Mach number contours In the 
flowfleld at the windward, waterline, and leeward rays are presented In 
figure 20. On the windward ray, the cone shock fall Inside the cowl lip, while 
on the leeward ray, the cone shock falls forward of the lip. The static pres- 
sure again Indicates the proper shock translation with the leeward shock mov- 
ing forward and becommlng stronger than the windward Interaction. Also evident 
are the expansions near the throat and Its Interaction with the shock system. 
The Mach number contours Indicate this same behavior, but also shows the bound- 
ary layer variations through the calculation. Work Is continuing on these 
three-dimensional shock wave boundary layer Interactions. 


Three-Dimensional Elliptic Invlscld/Vlscous Flow Solver (MINT) 

Specification of high grid resolution In the vicinity of a no-slip wall 
Is obviously required to resolve the wall boundary layers. However, specifi- 
cation of a high resolution region to resolve a normal shock wave Is not a 
simple problem because the shock location and shape may not be known a priori. 
Thus, a viable shock tracking adaptive mesh strategy was developed to properly 
resolve the regions of high gradients for normal shock wave turbulent boundary 
layer Interactions. In particular, a search for the maximum pressure gradient 
location was used to establish a definition for the shock wave center. Since 
any spurious oscillations In the solution could prove detrimental to accurate 
shock center determination, a “filtering scheme" was applied to the process of 
the shock location. Essentially, the procedure Identifies turning points In 
the wall static pressure distribution, examines the change In pressure between 
subsequent turning points and carries out a search for the maximum pressure 
gradient In the Interval having the largest pressure rise. Once the shock 
center Is located, a new grid Is constructed by centering a slnh function at 
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the new shock center location. The grid motion Is accounted for In the govern- 
ing equations through Inclusion of terms containing derivatives of the compu- 
tational coordinates with time. 

The present effort to benchmark verify the HINT general Invlscld/vlscous 
three-dimensional flow solver (fig. 21) has centered around the normal shock 
turbulent boundary layer Interaction In the channel of Mateer and Vlegas (ref. 
38) and In the variable area diffuser of Sajben, Bogar, and Kroutll (ref. 41). 
in particular, two main Issues were under Investigation, namely, (1) does an 
adaptive mesh or grid clustering region which moves with the shock region dur- 
Ing the course of the calculation Introduce any significant errors, and (2) 
what role does artificial dissipation play In solution accuracy. 

C onstant area two-dimensional normal shock wave Interaction . - The first 
case considered simulates the normal shock wave turbulent boundary layer In- 
teraction which will occure at modest upstream supersonic Mach numbers In a 
constant area tube with circular cross section. The case has been specified 
to match the experimental data of Mateer and Vlegas (ref. 13). The calcula- 
tion was performed for an Inlet Mach number of 1.44, an Imposed exit to Inlet 
static pressure ratio of 2.0, and a Reynolds number of 5.83E5 based on up- 
stream boundary layer thickness of 2.5 centimeters. In order to maintain ade- 
quate mesh resolution In all parts of the flow field, the adaptive mesh 
strategy described earlier was Invoked. The calculation was performed on a 
computational mesh with 41 transverse and 31 streamwlse points, using a mixing 
length turbulence model. In the original analysis paper by Roscoe, Shamroth, 
Glbellng, and McDonald (ref. 13) comparison between experimental data and 
analysis were made with both a mixing length and k-e turbulence model. The 
results obtained with the mixing length turbulence are shown In figures 22 and 
23. Figure 22 presents a comparison between the measured and calculated static 
pressure distribution while figure 23 presents the comparlslons for the stream- 
wise velocity profiles through the Interaction region. It Is apparent that 
the mean flow quantities are well predicted by this method and the adaptive 
mesh strategy described Introduced no appreciable error In the calculation. 

Variable area two-dimensional normal shock wave Interaction . - Figures 
24(a) and (b) present a comparlslon between MINT Navler Stokes calculations of 
the normal shock wave turbulent boundary layer Interaction In a variable area 
diffuser with the experimental data of Sajben, Bogar, and Kroutll. Calcula- 
tions were performed on a computational mesh consisting of 25 streamwlse points 
and 25 transvers nodal points. The computations were performed on the Lewis 
IBM 370/3033 high speed computer, and required 0.013 CPU secs per grid point 
^ e '' ^ me step. Since that calculation, the computing time has been reduced to 
0.0065 CPU secs per grid point per time step as a result of code "clean-up". 

The equlvelant time on the Lewis CRAY I high speed computer would therefore be 
6.5E-4 CPU secs per grid point per time step after vectorl zatlon. The calcu- 
lated top wall distribution Is shown In figure 24(a) and the calculated bottom 
wall distribution Is shown In figure 24(b). Solutions were obtained with two 
values of artificial dissipation, a = 0.05 and a = 0.5. It Is obvious that 
the choice of the artlflcal dissipation parameter significantly affects the 
calculated results. The normal shock Is captured properly with o = 0.05, 
while the results were severlly smeared with a = 0.5. In fact, the a =0.5 
calculations did not even contain a supersonic region (fig. 24). 

The excellent results obtained using the MINT three-dimensional flow 
solver established clearly that an adaptive grid approach can be used to 
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calculate shocked flow fields in which the position of the shock wave Is un- 
known a priori. Also, the role of the artificial dissipation parameter was 
examined In this study. For an Inappropriate choice of artlflcal dissipation, 
the solution will be severlly smeared and may, In fact, even suppress the 
appearance of the terminal shock region. A detailed discussion concerning the 
manner In which these calculations were performed as well as a more In depth 
examination of the effects of artlflcal dissipation appear Ip two publications 
by Liu, Shamroth, and McDonald (refs. 14 and 15). Future Investigation on the 
normal shock wave turbulent boundary layer Interaction will center on the 
problems consernlng the responce of the terminal shock to an externally applied 
disturbance and the effects of turbulence modeling on the small scale flow 
properties. Early efforts to understand the dynamic effects of a back pressure 
disturbance on the normal shock wave turbulent boundary layer Interaction 
appears In a paper by Liu, Shamroth, and McDonald (ref. 16). 


Three-Dimensional Subsonic Viscous Marching Flow Solver (PEPSIG) 

Recent effort to benchmark the PEPSIG general three-dimensional subsonic 
viscous marching flow solver has centered around the development of secondary 
flow In a number of subsonic duct configurations (fig. 25). The philosophy 
underlying this list of benchmark phenomena was to verify the PEPSIG three- 
dimensional flow field solver on a sequence of flows with Increasing complex- 
ity. These Include laminar and turbulent flow In a 90° bend with a square 
cross section of Taylor, Whltelaw, and Ylanneskls (ref. 44) laminar and tur- 
bulent flow In the 90° bend with circular cross section of Enayet, Gibson, and 
Taylor (ref. 45) laminar flow In the 180° circular pipe of Agrawal, Talbot, 
and Gong (ref. 46) laminar and turbulent flow In the S-shaped ducts with square 
and circular cross sections of Taylor, Whltelaw, and Ylanneskls (refs. 47 and 
48) the turbulent flow In the 45-45 circular S-shaped duct of Bansod and 
Bradshaw (ref. 49) and the turbulent flow In the square-to-round transition 
duct of Taylor, Whltelaw, and Ylanneskls (ref. 50). In addition, NASA Lewis 
also Is sponsoring an experimental effort at the University of Tennessee Space 
Institute to study the structure of secondary flows In S-ducts with high 
entrance Mach number and diffusion. Experimental results of this effort appear 
In a paper by Vaklll, Wu, Bhat, Liver, Hlngst, and Towne (ref. 51). 

Circular 90° Bend . - Extensive calculations were made for the flow geom- 
etry In which detailed measurements were made by Enayet, Gibson, and Taylor 
(ref. 45). This geometry consisted of a circular duct with a 90° circular- 
arc-bend and with straight sections both upstream and downstream of the bend. 
The ratio of bend radius to duct width was 2.3. The measurements were taken 
for Reynolds number of 790 (laminar flow) and 40 000 (turbulent flow). Two 
computational grid systems were used for the laminar evaluation; a coarse mesh 
system composed of 20X20X75 (30 000) nodal points and a fine grid system com- 
posed of 40X40X75 (120 000) points. On the Lewis CRAY I high speed computer, 
these two cases used 0.7 and 2.9 min of CPU time. Likewise, two grid systems 
were also used to evaluate the turbulent data; a coarse mesh composed of 
25X25X86 (53 750) points and a fine mesh system having 50X50X86 (215 000) nodal 
points. Calculations were performed on the Lewis CRAY I computer using a total 
of 1.3 and 5.1 min of CPU time respectively. A comparison between the experi- 
mental data and the analysis using the coarse and fine grid system Is pre- 
sented In figure 26 for the laminar flow case. Two regions of separation were 
encountered for the laminar flow case. The first separation region was 
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located at the entrance to the bend on the outer wall and the second on the 
Inner wall near the bend exit. A comparison of the measured and computed 
streamwlse velocity profiles (fig. 26) show that the "separation model" within 
the PEPSIG flow solver can simulate the effects of "weak" separations on the 
main flow field. Figure 27 presents the structure of the secondary flow at 
the four experimental measuring stations. The secondary flow Is formed because 
the fluid near the flow axis, having higher velocity, Is acted upon by a larger 
centrifugal force than the slower fluid near the walls. The faster moving 
fluid at the center moves outwards, pushing the fluid In the boundary layer at 
the outer wall around towards the Inner wall. Thus fresh fluid Is continually 
being brought Into the neighborhood of the outer wall and forced towards the 
Inner wall. The overturning that occurs In the boundary layer generates a 
strong vortex system, which migrates away from the wall near the bend exit. 
Presented In figure 28 Is a comparison between the measured and computed turb- 
ulent streamwlse velocity profiles at the four measuring stations through the 
90° bend. For turbulent flow, the results were very sensitive to mesh resolu- 
tion In the regions of high shear, as can be seen In figure 28. Figure 29 
presents the turbulent secondary flow structure which consists of a pair of 
counter-rotating vortlcles formed from the overturning of the flow within the 
wall boundary layer. This "overturning" can also be seen In the surface oil 
film patterns presented In figure 29. 

Circular 180° Bend . - Agrawal, Talbot, and Gong (ref. 46) have obtained 
detailed LDV measurements of laminar flow development In curved pipes with 
uniform entry velocity. This experimental data was used by Towne (ref. 53) to 
evaluate and verify the ability of the PEPSIG three-dimensional viscous flow 
solver to quantitatively predict the generation of pressure driven secondary 
flows in curved ducts. Experimental data was obtained at a Reynolds number of 
1263, based on the cross sectional radius, and an entry velocity corresponding 
to a Dean number of 565. A computational mesh consisting of 50X50 nodal points 
In the transverse plane and 226 forward marching steps, for a total of 565 000 
mesh points, was used for the results shown in figure 30. On the Lewis CRAY I 
computer, this calculation took 17 min of CPU time. A comparison between the 
measured and calculated streamwlse contour plots presented In figure 30 demon- 
strate that for well prescribed geometry description and initial data, the 
PEPSIG flow solver can simulate the fine detail of flow structure associated 
with developing pressure driven flow fields. 

Circular 22.5-22.5° S-Bend . - Taylor, Whltelaw and Ylanneskls (ref. 48) 

In an experimental Investigation sponsored by NASA Lewis Research Center, ob- 
tained a series of three-component LDV velocity measurements on the structure 
of the flow that develops In a 22.5-22.5° circular S-bend. Measurements Include 
laminar flow at a Reynolds number of 790 and turbulent flow at a Reynolds 
number of 48 000. Two computational grid systems were used for the laminar 
evaluation; a coarse mesh system composed of 20X20X80 (32 000) nodal points 
and a fine grid system composed of 40X40X80 (128 000) points. On the Lewis 
CRAY I high speed computer, these two cases used 0.8 and 3.0 min of CPU time. 
Likewise, two grid systems were also used to evaluate the turbulent data; a 
coarse mesh composed of 25X25X80 (50 000) points and a fine mesh system having 
50X50X80 (200 000) nodal points. Calculations were performed on the Lewis CRAY 
I computer using a total of 1.2 and and 4.8 min of CPU time respectively. 

Figure 31 presents a comparison between the measured and computed laminar 
development of streamwlse velocity profiles at the four measurement stations 
through the S-duct. Because the boundary layers with this duct were large 
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relative to the duct dimensions, the laminar flow calculations were not sensi- 
tive to mesh resolution. The complex secondary flow structure that can develop 
In S-duct configurations Is graphically shown In figure 32. The system of 
counter rotating vortex pairs develops In the first bend, and In the second 
bend, the secondary flow begins to reverse forming another pair of counter 
rotating vortlcles. Also presented In figure 31 are the wall shear stress 
signatures that result from this secondary flow structure. Much greater sen- 
sitivity to resolution of high shear regions was encountered for the turbulent 
flow development through the S-duct (fig. 33). Qualitatively, the turbulent 
flow development resembles the laminar field; however, the high Reynolds number 
In the developing turbulent flow results In less severe secondary flow (fig. 

31). Again, the computed results are In very good agreement with both the 
laminar and turbulent measurements (ref. 47). 

Circular 45-45° S-Bend . - A series of measurements obtained at the 
Department of Aeronautics, Imperial College of Science and Technology, by 
Bansod and Bradshaw (ref 49) were used to verify the PEPSIG flow solver at 
higher Reynolds numbers. The configuration used for this study was the 45-45 
symmetric short Intake S-shaped duct with an R/D of 2.25. The duct entry 
velocity was nominally set at 45 meters/sec, which gave a Reynolds number of 
5.0E5 based on duct diameter. Measurements were presented of total pressure, 
static pressure, surface shear stress, and yaw angle for the flow through the 
S-shaped duct. The computational mesh used for this study was 50X50X100, which 
gave a total number of grid points of 250 000. Computations were performed on 
the Lewis CRAY I high speed computer using 6.0 min of CPU time. Shown In 
figure 35 Is a comparison between the calculated and measured surface shear 
stress at three circumferential surface lengths labeled the N-length, E-length, 
and S-length. Excellent agreement was obtained In spite of the fact that a 
simple eddy viscosity turbulent model was used In this calculation. Also shown 
In figure 35 Is the surface shear stress color signature for this S-shaped 
duct. The small region of separation or near separation that was observed by 
Bansod and Bradshaw along the N-length Is clearly visible from the wall shear 
stress signature. Figure 36 shows a comparison between the calculated and 
measured total pressure loss contours at the compressor face station In addi- 
tion to the secondary velocity vector flow field. It Is clear that the PEPSIG 
three-dimensional flow solver captured the proper flow physics at the com- 
pressor face Including the pair of counter-rotating vortices In the boundary 
layer observed (ref. 49). 

Square-to-round transition duct . - Included In the series of experiments 
sponsered by NASA Lewis Research Center at Imperial College of Science and 
Technology, London, was flow In a square-to-round transition duct. These re- 
sults have been published by Taylor, Whltelaw and Ylanneskls (ref. 50) and 
Include three-component LDV measurements of the velocity field through the 
transition duct. The flow In this straight centerline transition duct was 
turbulent at a Reynolds number of 35 350. The computational mesh used for 
this calculation was 25X25 nodal points In the cross plane and 51 forward 
marching stations for a total number of 31 875 grid points. On the Lewis CRAY 
I computer, the CPU time was less than 1.0 min. Comparisons between the meas- 
ured and calculated streamwlse velocity profiles In a symmetry plane are pre- 
sented In figure 37 at each reported measurement station through the transition 
duct. This comparison shows very good agreement between analysis and 
experiment . 
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Excellent results have been obtained using the PEPSIG three-dimensional 
viscous marching flow solver to model the development of pressure driven sec- 
ondary flows, using a number of basic benchmark data sets for evaluation. 
Extensive evaluations were made for laminar flow within well defined duct 
geometries and well documented Initial data. The excellent agreement with the 
laminar data demonstrates the accuracy of the solution algorithm without the 
uncertainty of describing the turbulent transport process by means of a turbu- 
lence model. This demonstration of accuracy should be a prerequisite to de- 
termining differences between numerics and physics, and In particular, the 
effects of turbulence model on solution results. In general, the flows studied 
were dominated by pressure forces rather than shear forces, hence, the Influ- 
ence of normal stress driven flows was small, and the effective viscosity 
approach was appropriate. Even with "weak" separations, the flows were pri- 
marily pressure controlled, although turbulent mixing became more Important. 
Still, the effective viscosity approach may likely be the best compromise, but 
consideration should be given to stress turbulence models. These evaluations 
have established the Importance of mesh point resolution In calculating the 
structure of secondary flow In subsonic diffusers, and consequently diffuser 
performance, compressor face recovery, and distortion. This presupposes that 
the duct geometry Is well defined and computationally smooth, In addition to 
having an accurate description of the Initial data. 

Rockwell BIB Inlet duct . - Figures 38 and 39 present the results of a 
series of calculations on the Rockwell BIB Inlet duct using the PEPSIG three- 
dimensional viscous marching flow solver. The purpose of the computer study 
was twofold; namely (1) examine the difficulties associated with surface fit- 
ting a "real" Inlet duct geometry and (2) determine the effects of geometry 
description on the accuracy of the solution. The Rockwell BIB Inlet duct 
transitions from nearly rectangular cross-section at the Inlet face to a cir- 
cular cross-section at the compressor face, with a double S-bend Insert be- 
tween these two stations. The computations were Initiated with an entrance 
Mach number of 0.5 and Reynolds number of 1.0E6 per ft. Calculations were 
performed with a computational mesh consisting of 50X50X100 (250 000) nodal 
points using 7.5 min of CPU time on the Lewis CRAY I computer. Shown In 
figure 38 Is the computed static pressure signature on the Inside surface of 
the BIB Inlet duct along with the surface analytical oil film patterns. The 
surface oil film patterns highlight the over turning that results In the double 
S-bend Insert portion of this Inlet duct, which helps form the total pressure 
distortion presented In figure 39. The "quality" of this calculation depended 
very strongly on the ability to describe a computationally "smooth" Inlet duct 
surface and to compress the mesh Into regions of high shear. Complex surface 
geometries, such as the BIB Inlet duct, are very difficult to describe compu- 
tationally, and much time and effort had to be expended on this problem. How- 
ever, care In constructing the geometry and mesh system In the Initial steps 
of a duct analysis pays dividends later In the design-analysis process. 


USER ORIENTED CODE DEVELOPMENT AND PROTOCOL STUDIES 

The results of the comparisons between the PEPSI/MINT analysis and bench- 
mark experimental data demonstrate a strong need to establish a computation 
protocol to Insure the analytical results are reasonable and the desired In- 
formation correct. Careful attention has been paid to those factors which 
could substantially affect solution accuracy. In particular, the need for 
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adequate mesh resolution In regions of high shear has been shown to be a crit- 
ical factor affecting solution accuracy In many flow problems. Use of exces- 
sive artificial dissipation to suppress spatial oscillations In the solution 
has also been shown to suppress critical physics. Accurate and computationally 
"smooth" duct geometries along with well defined Initial data are also criti- 
cal factors affecting solution accuracy. Thus the results from benchmark 
validations of the PEPSI/HINT codes are being carefully considered to assess 
their meaning and Implications with particular emphasis on sensitive control- 
ling parameters. Based on these studies, formal "case running" protocols or 
structure for code application will be developed to guide the general user 
community. 

Another primary Item to be considered under this effort would focus upon 
the use of computer graphics In understanding the highly complex flows that 
develop within advanced supersonic Inlet systems. The use of graphics In these 
areas has been under development for several years at Lewis Research Center 
Including advanced techniques for visualizing complex three-dimensional flows 
In Inlet, nozzle, and mixer flow fields using color graphical techniques 
( refs . 54 and 55) . 


THREE-DIMENSIONAL GRID GENERATION AND NUMERICAL TECHNIQUES 

The computational mesh Is one of three key factors that have a profound 
effect on the accuracy of the numerical solution; the other two factors being 
the numerical solution algorithm and the flow field Itself. These factors 
cannot be separated. In each of the zonal regions within the Inlet system, 
the construction of the computational grid Is the key factor. However, the 
relevant physics to be studied In each of these propulsion components may re- 
quire elliptic solutions, either Euler or Navler Stokes, for all or part of 
the flow field. In turn, the construction of the computational mesh cannot be 
performed Independently of the available solution algorithm. There are no 
quantifiable criteria for generating acceptable finite-difference grids, 
Including an assessment of grid smoothness, skewness and cell aspect ratio 
Independently of the numerical solution obtained on the mesh. Lastly, If 
three-dimensional elliptic solutions are to become routinely used, even on 
advanced computers such as the CRAY, then major efforts are required to In- 
crease the computing speed and effectiveness of the solution algorithm. 

In solving the equations of fluid dynamics by numerical means, two major 
problem areas are encountered. The first Is three-dimensional mesh generation 
and the second Is the algorithm by which the governing equations are processed 
to obtain a solution. Mesh generation Involves three distinct processes, (1) 
the description of the bounding surfaces, (2) construction of the grid lattice 
within the bounding surfaces, and (3) the solution dependent adaptation or 
clustering of the grid lattice to accurately define the region of rapid solu- 
tion variation. The numerical algorithm used to obtain the solution also In- 
volves three distinct processes, (1) replacement of the system of equations 
with a discrete representation, (2) stable and efficient solution advancement, 
and (3) data base management to control the flow of Information to and from 
the numerical solution generator. 
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Three-Dimensional Mesh Generation 


Description of bounding surfaces . - Looking at each of the three processes 
Involved In mesh generation, It Is noted that a considerable amount of litera- 
ture has evolved on the fitting of surfaces for CAD/CAM and other applications. 
Much use of this existing literature could be made for fluid flow problems, 
but to date It has not been exploited. It Is recommended that some effort be 
expended to review, assess, and explore the possible use of this existing 
technology In the fluid mechanics of three-dimensional mesh generation. 

Construction of the grid lattice . - The construction of the grid lattice 
within and fitted to the bounding surface can be accomplished In several ways. 
Most notable, the grid lattice can be constructed numerically either as the 
solution to a partial differential equation or geometrically as a parameteri- 
zation of the bounding surfaces. Either way, the grid lattice construction 
requires some trial and error on the part of the analyst and the range of con- 
figurations Is often limited. Fundamental studies of both the aforementioned 
approaches to grid construction are required to understand and Improve the 
accuracy, reliability, generality, and ease with which grids may be generated. 
Of particular note Is the development of mesh smoothness criteria and param- 
eters which Influence this smoothness, leading ultimately to Improved grids 
for fluid dynamics computations. Most Importantly, while basic understanding 
can be obtained In one or two space dimensions, this Investigation will empha- 
size three-dimensional grid generation for the type of fluid flow problems 
encountered In Inlet, mixer, and nozzle components. 

Adaptive mesh clustering . - Adaptive mesh clustering Is a very powerful 
technique to effectively use a limited number of grid points to define the 
multiple length scales which can arise In fluid dynamic problems, such as 
boundary layer and Invlscld core regions. In numerous problems, particularly 
In three-dimensional space, the location of the region requiring definition Is 
not known a priori, but emerges as the solution develops. A transonic shock 
wave boundary layer Interaction would be a case In point. For such problems, 
adaptive mesh clustering Is a very powerful technique. Particularly In three 
dimensions, sufficient mesh resolution cannot be afforded without some mesh 
clustering and knowledge of the extreme solution variation. Solution depen- 
dent mesh adaptation can cause Instability and great care Is required to obtain 
stable, accurate, and, In some cases, optimal mesh clustering. Consequently, 
basic studies of the solution dependent mesh adaptation algorithms are required 
to understand and Improve the stability, accuracy, distribution, and adaption 
process In three-dimensional space. 


Numerical Algorithm Development 

Higher order discrete representation . - In the numerical solution of 
multidimensional flow problems, particularly In three dimensions, a major 
limitation Is the available computer memory. To the user this means that the 
total number of grid points In the computational space must be limited with a 
concomitant degradation In the solution accuracy. Conversely, when a desired 
level of accuracy Is sought, In many Instances the required number of spatial 
grid points needed to achieve that goal could be excessively large so as to 
make the calculation either Impractical or extremely expensive. Hence, there 
Is a great Interest In minimizing the total number of grid points and accruing 
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the associated benefits which Include reduction In computing time, a reduction 
In core requirements and a reduction of data transfer. All of these factors 
contribute to lower computing costs. 

A method of achieving these desired benefits Is the use of higher order 
spatial difference approximations. The replacement of the continuous system 
of equations by an accurate discrete representation Is especially Important In 
three-dimensional space.' In one dimensional space It Is known that a more 
accurate discrete representation of the governing equations can be obtained 
which results In a reduction In the number of grid lattice points required for 
an accurate solution. In multidimensional space, the required grid point 
density In one dimension Is raised to the power of the space dimension. This 
raising to the power transforms the modest grid point savings In one dimen- 
sional space to Impressive benefits In three-dimensional space. Thus more 
accurate higher order discretization schemes should be Investigated from the 
view point of multidimensional use. 

Solution advancement . - The second process, that of solution advancement, 
has two major problem categories; (1) time dependent or unsteady flows and (2) 
steady state flow problems. Fluid flow problems that require the transients 
to be resolved, such as subcrltlcal self-sustained oscillations (Inlet buzz), 
Introduce considerations which are not present If only a steady state solution 
Is sought. The transients generally pose a hyperbolic problem In time which 
Is further complicated If the flow Is Incompressible or at a very low Mach 
number. In addressing the problem of accurately and efficiently resolving the 
transients, the needs and limitations of three space dimensions must be kept 
In mind. Basic studies would be expected to lead to more robust and economical 
time dependent multidimensional solution algorithm for subsonic, transonic, 
and supersonic flows. 

The steady state flow problem of Interest poses an "elliptic" or "mixed" 
problem In which all the boundary conditions affect all or nearly all of the 
Interior flow field. The problem Is also nonlinear, thus some form of Itera- 
tion may be required. This In many schemes Is added to a basically Iterative 
approach to solving the elliptic problem. This has led to the Investigation 
of the "fast" Iterative elliptic solvers for some problems and fast direct 
schemes for other more restrictive problems. The fast "Iterative" schemes are 
more general and are applicable to the Navler Stokes equations, the "parabo- 
lized" Navler Stokes equations and the Euler equations and hence should be 
actively pursued. A number of schemes have shown considerable promise on 
single equations such as matrix pre-conditioning, relaxation parameter cycling 
and even mesh cycling (multi-grid techniques) and have demonstrated Improve- 
ment In the rate of convergence to steady state solution. The role of con- 
vergence acceleration In systems of equations Is presently not well understood 
and should be actively Investigated. 

Data base management . - The third consideration for the solution algorithm 
Is the data base management aspects of the problem. Here the architecture of 
the computer Itself enters from the point of view of the available memory, 
data transfer capability, and processor parallelism. The CRAY computer Is 
fast becoming the Industry standard for large scientific computing and conse- 
quently It Is advocated that this type of computer be used to guide the data 
base management aspects of the problem. Vectorl zatlon of the solution algo- 
rithm should be examined, but also In view of the large arrays Involved In 
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three-dimensional flow calculations, the algorithm data and data transfer needs 
should be considered. It Is possible that a highly structured algorithm with 
local data requirements could out perform another algorithm on the CRAY simply 
by vlrture of modest storage data transfer requirements. Thus, the considera- 
tion of the solution algorithm and the computer architecture must be borne In 
mind when evaluating the relative efficiency of the various schemes. 


CONCLUDING REMARKS 

A general system of three-dimensional viscous marching and elliptic 
solvers called PEPSI/MINT has been developed at the NASA Lewis Research Center 
for the analysis and design of advanced technology supersonic Inlet systems. 
Through a systematic and careful evaluation of these solvers by a two tier 
validation process Involving generic verification and detailed benchmark vali- 
dation data, the applicability and performance of these codes Is being assessed 
and documented. From the benchmark validation process, a formalized protocol 
Is being developed to assist the user In applying these three-dimensional flow 
solvers to Insure reasonable correctness and to be able to distinguish between 
numerics and physics. In addition, an extensive library of plotting and 
three-dimensional graphical routines were developed to assist the user In 
understanding the very complex flow fields emerging as solutions from these 
solvers . 

Major advances In computational fluid dynamics have emerged from this 
program, but significant problems remain and these are being addressed. In 
particular, future effort will be directed towards the development of Improved 
three-dimensional grid generation, numerical methodology, and computer tech- 
nology utilization. Likewise, the detailed benchmark validation of the 
PEPSI/MINT solvers will continue considering more complex flow fields. 
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APPENDIX - GENERALIZED THREE-DIMENSION FLOW SOLVERS 

An extensive effort Is currently underway to structure the PEPSI/MINT 
three-dimensional flow solvers described In this section, such that they have 
the ability to communicate with each other. This Is In keeping with the 
philosophy of the zonal methodology concept. Common file storage for all 
dependent variables, Identical formatted restart Input/output files, uniform 
geometry methods for all- PEPSI/MINT codes and common plot files, to Interface 
with the extensive library of existing three-dimensional plotting and color 
graphical programs (refs. 54 and 55) are among the tasks being undertaken. 


Three-Dimensional Supersonic Viscous Marching Flow Solver (PEPSIS) 

The PEPSIS code that has been developed (refs. 1 to 5) solves the set of 
three-dimensional forward marching equations applicable to flow fields within 
supersonic Inlets and nozzles of axlsymmetrlc, two-dimensional or conformal 
geometries. The equations used In the analysis are based on the parabolized 
form of the three-dimensional, steady state, ensemble-averaged Navler Stokes 
equations. In the supersonic region of the core flow field, the required 
parabolizing assumption simply reduces to the neglect of the streamwlse diffu- 
sion terms. In the flows anticipated for this region, streamwlse diffusion Is 
negligible and, therefore, the assumption Is not restrictive. In the subsonic 
portion of the flow field, further assumptions are necessary to control the 
appearance of branching solutions. Suppression of branching and development 
of a stable forward marching procedure requires replacing the normal momentum 
equation In the subsonic region by the usual boundary layer approximation, 
which, for zerc curvature, sets the normal derivative to zero and for curved 
walls, allows for coordinate curvature effects. In addition, the wall tangency 
condition Is applied to the entire subsonic region. This latter condition 
replaces the continuity condition at the sonic line by a specified flow direc- 
tion. Since the subsonic region for most problems of practical Interests Is 
expected to be small, these approximations are not expected to have detremental 
effects In calculating Inlet flow fields. 

The PEPSIS three-dimensional flow field solver has been applied to a 
variety of two and three-dimensional flow problems at supersonic speeds In the 
absence of a major separation zone. The flow solver accepts a general ortho- 
gonal coordinate system through an external mesh generator which generates 
computational grid points In the given physical domain. Boundary conditions 
are set via Input data flags. Among those available are function, first de- 
rivative, and second derivative of the velocity, density, and temperature, 
specification of static pressure with Its derivatives, momentum equation, 
one-sided Invlscld governing equations, Mach line extrapolation, and wall 
function formulation. Inclusion of the energy equation as well as Its form 
and the choice of two or three dimensions and viscosity options are controlled 
via Input data flags. The PEPSIS flow field solver has been validated through 
detailed comparison with a large variety of two and three-dimensional experi- 
mental flow fields (refs. 22 to 37). 


Three-Dimensional Elliptic Invlslcd/Vlscious Flow Solver (MINI) 

The approach being used within the Invlscld/vlscous MINT flow field solver 
applies the consistently split Linear Block Implicit (LBI) procedure (refs. 6 
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to 16) to the Navler Stokes equations. The method can be briefly outlined as 
follows: the governing equations are replaced by an Implicit time difference 
approximation, optionally a backward difference or Crank-NIcolson scheme. 

Terms Involving nonllnearltles at the Implicit time level are linearized by a 
Taylor series expansion In time about the solution at the known time level, 
and spatial difference approximations are Introduced. The result Is a system 
of multidimensional coupled, but linear, difference equations for the dependent 
variables at the unknown or Implicit time level. To solve these difference 
equations, the Douglas-Gunn procedure (ref. 10) for generating Alternating 
Direction Implicit (ADI) schemes as perturbations of fundamental Implicit dif- 
ference schemes Is Introduced. This technique leads to systems of coupled 
linear difference equations having narrow banded matrix structure which can be 
solved efficiently by a standard block elimination method. 

The HINT three-dimensional flow field elliptic solver accepts a very 
general coordinate system through Input of the physical location of the compu- 
tational grid points. Boundary conditions are set through Input data flags. 
Among those available are function, first derivative and second derivative of 
the velocity, density and temperature, specification of static and total pres- 
sure or their derivatives, and one-sided governing equations. Inclusion of an 
energy equation as well as the choice of two or three dimensions and turbu- 
lence model are controlled via Input data flags. Specification of very high 
Reynolds number along with Implementation of proper boundary conditions would 
allow the elliptic solver to be applied to the Euler equations. Finally, out- 
put Including flow field arrays and generation of plot files are also con- 
trolled via Input data flags. The HINT flow solver has been validated through 
two and three-dimensional calculation of transonic diffuser flows with a 
terminal shock wave (refs. 33 to 37). 


Three-Dimensional Subsonic Viscous Harchlng Flow Solver (PEPSIG) 

The PEPSIG three-dimensional flow solver (refs. 17 to 21) Is a general 
approach used for predicting subsonic flows In three-dimensional passages hav- 
ing little or no streamwlse separation and Is based on the primary-secondary 
velocity decomposition method of Briley and HcDonald (ref. 16) for application 
to viscous subsonic flow In smoothly curved geometries. The objective of this 
approach Is to Introduce approximations which adequately represent essential 
physical features of Interest and yet lead to governing equations which can be 
solved much more economically than Navler Stokes equations. In the present 
application. It Is necessary to provide an adequate representation of primary 
flows, secondary flows, viscous effects, and their local Interactions. An 
Invlscld flow solution Is first obtained for the geometry In question. The 
Invlscld flow satisfies an elliptic governing equation requiring downstream 
boundary conditions and thus Includes transverse variations In streamwlse 
pressure gradients usually associated with flows In curved passages. This a 
priori pressure field Is then used as an Imposed criteria upon a set of equa- 
tions which constitutes a well-posed Initial-value problem In space. Solution 
of these equations leads to a prediction of both the velocity field and the 
corrected pressure field. When solved as an Initial value problem, the tech- 
nique provides a reduction In computational effort of one or more orders of 
magnitude over Navler Stokes equations. 
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The PEPSIG three-dimensional flow field solver Is applicable to complex 
geometries having curved and twisted centerlines, variable cross-sectional 
areas, and shapes which require the use of nonorthogonal body fitted coordi- 
nate systems. Quantitative assessment of the method's predictions has been 
made by comparison with experimental data for a variety of geometries and 
fluid dynamic conditions (refs. 43 to 53). 
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ANALYSIS 
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Figure 1. - Supersonic inlet zonal methodology. 
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Figure 3. -PEPSIS benchmark verification. 
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Figure 4 - Two dimensional oblique shockwave turbulent 
boundary layer interaction; static pressure distribution. 
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Figure 5. - Two dimensional oblique shock wave turbulent bound 
ary layer interaction, Mach number profiles. 
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Figure 6. - Three dimensional sidewall shock wave turbulent boundary 
layer interaction: yaw angle distributing. 
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Figure 9. - Concluded. 
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Figure 11. - PEPSIS benchmark verification, Lewis 1X1 foot supersonic 
wind tunnel experiments. 




Figure 12. -Three dimensional glancing sidewall/corner shockwave interaction, analytical surface oil film patterns. 
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gure 13. - Three dimensional glancing sidewall/corner interaction 
experimental surface oil film patterns. 


PHENOMENON 


TYPE OF DATA 


SOURCE 


DETAIL BL SURVEYS 


CRESCI 


• HYPERSONIC CORNER 


• INTERSECTING WEDGE WEST AND 

KORKEGI 


DETAIL BL SURVEYS 


DETAIL BL SURVEYS 


• CROSSED SIDEWALL 
SHOCK/BL INTERACTION 


PROPOSED 


• SKEWED WEDGE 


SETTLES AND DETAIL BL SURVEYS 
BOGDONOFF 


Figure 14 - PEPSIS benchmark verification for advanced air 
craft inlet concepts. 



Figure 16. - Hypersonic corner interaction, Mach number signature. 






Figure 17. - Hypersonic corner interaction, flow field total pressure 
signature. 
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Figure 19. - Lewis 40/60 mixed compression supersonic inlet, surface 
static pressure, Moo=3.0, a = 2.0°. 
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Figure 20. - Lewis 40/60 mixed compression supersonic inlet, flow field 
static pressure and mach number Me =3.0, a = 2.0°. 
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Figure 21. - MINT benchmark verification. 
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Figure 22. - Constant area normal shock wave interaction, static 
pressure. 
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Figure 23. - Constant area normal shock wave interaction velocity 
profiles. 
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(b) Top wall static pressure. 


Figure 24. - Variable area normal shockwave interaction. 
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Figure 25. - PEPSIG benchmark verification, 
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Figure 26. - Circular 90° bend, laminar flow, streamwise velxity. 
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Figure 29. - Circular 90° bend, turbulent flow, secondary velocity. 




STREAMWISE VELOCITY 


Figure 30 . - Circular 180 ° bend, laminar flow. 



Figure 31. - Circular 22.5-22.5° S -bend, laminar flow, streamwise 
velxity. 




Figure 33. - Circular 22. 5-22. 5° S-bend turbulent flow, streamwise 
velocity. 
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Figure 35. -Circular 45-45° S-bend, turbulent flow, wall shear stress. 
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Figure 37. - Square-to-round transition duct, turbulent flow, 
streamwise velocity. 
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Figure 38. - Rockwell BIB inlet duct surface conditions. 
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